On continuous variable quantum algorithms for oracle identification problems 
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I. INTRODUCTION 

Quantum information protocols have been demon- 
strated experimentally in both the discrete- variable (DV) 
and so-called continuous-variable (CV) settings. DV 
quantum information protocols employ qubits 1| and qu- 
dits 0, and CV quantum information protocols regard 
continuously parameterized canonical position states as 
the logical elements analogous to qubits for the DV 
case [3(. CV quantum information is experimentally ap- 
pealing because sophisticated squeezed light experiments 
have led to claims of successful quantum information pro- 
tocols such as teleportation [J], key distribution [5j, and 
memory [(| 0] , but the theoretical status of CV quantum 
information is challenged by unresolved issues concerning 
quantum error correction [8(, non-distillability [9(, no-go 
theorems for quantum computation [ll[ , and the ab- 
sence of full security proofs for key distribution. 

CV information processing has also been studied for 
classical models L including the now named Blum-Shub- 
Smale machine fl2| and continuous Turing machines [l3| . 
These models are of background relevance to the research 
into CV quantum information and are referenced here for 
contextual purposes. 

In this paper, we establish a sound theoretical frame- 
work for studying quantum algorithms and apply this 
framework to study the CV analogue of the early DV 
quantum algorithm, known as the Deutsch-Jozsa (DJ) 
algorithm [Tj, [jja, • The problem solved by D J algo- 
rithm is the following. 

Problem 1. Given a function f : {0, 1}™ — > {0, 1} that 
is promised to be either constant (f takes the same value 
everywhere) or balanced (f takes the value on exactly 
half the inputs), determine whether f is constant or bal- 
anced. 

The best classical algorithm requires 2 n_1 + 1 evalu- 
ations in the worst case. If error is tolerated, for any 



integer m > 2, to achieve an error of at most 2~ m , any 
probabilistic algorithm requires a number of evaluations 
that is at least of order m [l7| • If the function is acces- 
sible on a quantum computer as a quantum oracle, then 
the DJ algorithm is exact and requires just one evalua- 
tion to solve the problem. 

Our focus here is on the CV analogue of the DJ algo- 
rithm, and we are inspired by the Braunstein and Pati 
formulation [l8[ of the CV DJ algorithm; however, our 
work differs from theirs in that ours relies only on logical 
states that are elements in the Hilbert Space and thus 
provides a strict CV version of the DJ problem. We in- 
troduce a particular model for the computation of the D J 
problem in a CV setting. Within the constraints of this 
model, our analysis shows that the CV DJ algorithm is 
necessarily probabilistic and its performance must there- 
fore be compared to the classical case where bounded 
error is tolerated and not to the classical deterministic 
case. 

We choose the DJ algorithm for the following reasons. 
Two types of quantum algorithms dominate the field, 
those that implement a version of the hidden subgroup 
problem and those that use a version of Grover's search 
algorithm [lLGil. An early example of the former is the 
Deutsch-Jozsa (DJ) algorithm [lj|, which is amongst the 
oracle class of problems [2Cj that have been important in 
demonstrations of quantum speed-ups. Finally the CV 
D J algorithm has a head start in the work of Braunstein 
and Pati so our analysis can build on their concepts [l|| . 

Our paper is presented as follows. In Sec. [Ill we review 
the DJ algorithm. Although this algorithm is well known, 
our review serves as a foundation for careful construction 
of the CV version of this algorithm. Furthermore we 
compare the DJ algorithm's performance against both 
deterministic and probabilistic strategies, especially be- 
cause the CV case can only be properly compared against 
probabilistic strategies because the CV DJ algorithm can 
never be deterministic. Our description of the DV D J al- 
gorithm comprises three steps so that these steps can be 
discussed separately during the construction of the CV 
analogue. 
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Our approach emphasizes a recasting of the DV D J al- 
gorithm in that we do not need the target qubit. This 
approach leads to an easier adaptation to the CV case. In 
Sec. UH we review the formalism of rigged Hilbert spaces 
(RHS) [2l[ since our CV algorithm, as well as any other 
CV quantum algorithms, must work in a RHS. This will 
have implications when we discuss the limitations of er- 
ror inherent in our CV DJ algorithm in Sec. IIIII and in 
Sec. ED 

In Sec. IIIII we adapt the DJ problem to the CV case 
and develop the CV DJ algorithm through the same three 
fundamental steps of the algorithm. We pay particular 
attention to the challenge of encoding a finite iV-bit string 
into functions over the real numbers. Overcoming this 
challenge enables us to recognize that perfect encoding 
results in the inability to determine if the encoding is of 
a constant string or balanced in a single execution of the 
algorithm. We show that this probabilistic nature of the 
algorithm is the result of an uncertainty relation between 
the continuous representation and its Fourier-transform 
dual representation. 

In Sec. IIV1 we determine an upper bound on the query 
complexity of the CV DJ algorithm. We note that be- 
cause the CV DJ algorithm is shown to be probabilistic, 
its performance can only logically be compared to the 
classical probabilistic algorithm and not to the classical 
deterministic algorithm. We conclude that the formalism 
presented herein is applicable to a wide range of oracle 
identification problems in a CV setting. 



with the promise that either f € So or f G Si, determine 
the index b such that f G Sf,. 

For N = 2™, we impose lexicographic order on the N- 
bit strings of {0, 1}™. We can then specify any function 
f z by writing all its N function values in a list z € {0, 1}^ 
of length N. The i th bit Zi in the list is 1 if / takes the 
value 1 on the i th bit-string of {0, 1}™. There are 2 N 
functions from n bits to one bit, and thus our candidate 
set has cardinality upper bounded by M < 2 N . In the 
following, we often write f z to denote the function that 
corresponds to the TV- bit string z. 

We are interested in finding an efficient strategy to 
identify the property of whether / belongs to set So or 
to set Si without necessarily determining / itself. In the 
DJ case, the property we are interested in is whether 
/ is balanced or constant [3, HE Gj| ■ The cost of the 
algorithm is the number of queries made to the oracle. 

With the promise of balanced or constant functions, 
there are far fewer than 2 N functions. The number of 
balanced and the number of constant functions is readily 
ascertained from the binomial theorem applied to power 
sets. The strings z of length N that correspond to the 
constant functions are the string consisting only of Is 
and the string consisting only of 0s. There are thus just 
two constant functions. The strings z of length N that 
correspond to balanced functions are the strings in which 
exactly half of the bits are and half are 1. There are 
thus precisely (^ 2 ) balanced functions. 



II. BACKGROUND 

We cast the DJ problem into the class of 'oracle iden- 
tification problems' in Subscc. Ill Al We then review de- 
terministic algorithms in Subsec. Ill Bl and probabilistic 
algorithms in Subsec. Ill CI In Subsec. Ill Dl we analyze 
an alternative representation of the quantum DJ algo- 
rithm that uses n qubits instead of the traditional n + 1 
qubits. In Subsec. [II El we present a primer on the rigged 
Hilbert space and close with a discussion of the concepts 
required to transition from discrete variables to continu- 
ous variables. 



A. The Oracle Identification Problem 

The DJ problem is an identification problem in which 
we are given a function from some candidate set S = 
{/li Hi ■ • ■ > /m} of functions. The candidate set S = 
SqUSi is the disjoint union of two collections of functions, 
and our task is to determine which of the two collections 
the function / is drawn from. 

Problem 2. Consider the set of all functions from n bits 
to one bit, T = {/ | / : {0, 1}™ -> {0, 1} }. Let S a and Si 
be disjoint subsets of T. Given some oracle 

/:{Q,ir-{0,l} (2.1) 



B. The Classical Deterministic Approach 

On a classical Turing machine, Problem[T]can be solved 
deterministically. A deterministic algorithm corresponds 
to submitting queries in the form of n-bit inputs and 
obtain the one-bit output for each query. There are N 
unique input strings, but the promise of balanced versus 
constant functions implies that only N/ 2 + 1 are required 
to determine whether the given function is balanced or 
constant, with certainty. 

The reason that fewer than N/2 + 1 queries is insuffi- 
cient is that only N/2 queries may reveal all output bits 
being the same, suggesting a constant function, whereas 
the remaining N/2 outputs could all be the opposite of 
the first N/2 queries. 



C. The Classical Probabilistic Approach 

In Subsec. Hi Bl we saw that fewer than N/2 + 1 queries 
is insufficient for a deterministic algorithm, but that case 
seems highly unlikely. More formally, fewer than N/ 2 + 1 
queries will identify most of the balanced functions as 
non-constant in much fewer than N/2 + 1 queries. Here 
we ask the question about how many queries are required 
if we are prepared to tolerate a small number of errors. 
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FIG. 1: Quantum circuit implementing the Discrete Variable 
DJ Algorithm. The upper line represents the n-qubit "con- 
trol" state, and the lower line represents the 1-qubit "target" 
state. 



In fact a probabilistic algorithm achieves an exponen- 
tially small error of 2~ m with a number of queries that 
is only linear in m [17| . To understand how a probabilis- 
tic algorithm can help, consider that, although a single 
query with a random input provides no information, two 
queries with two random inputs can be highly informa- 
tive. If the output from the second query differs from the 
first output, then the function is proved not to be con- 
stant and therefore must be balanced. If, on the other 
hand, the second output is the same as the first, then the 
outcome is not certain, but the more times the outputs 
are the same, the more confident one can be about the 
function being constant. 

We calculate the probability of successfully determin- 
ing whether the given function f z is balanced or con- 
stant. A lower bound on the success probability Pre- 
fer m queries can be achieved by examining a sampling- 
without-replacement strategy, which is expressed as 



p r =i n 



N/2 - (j - 1) 
N-(j- 1) 



> 1 



(2.2) 



Here the equality is calculated assuming sampling with- 
out replacement and shows dependency on N, whereas 
the inequality in Eq. (|2.2[) is based on sampling with re- 
placement and is independent of N. The failure proba- 
bility 1 — Piy declines exponentially in m, the number 
of queries. 

In Subsec. Ill Dt we study the quantum DJ algorithm 
next where we show that the problem can be solved with 
a single query independent of N. Although this exponen- 
tial speed-up is impressive when compared to the classical 
deterministic approach, it is much less so when compared 
to the classical probabilistic approach. 



D. The Quantum DJ Algorithm 

The quantum DJ algorithm has been shown to solve 
Problem [1] in a single query [l4T |. The quantum DJ algo- 
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FIG. 2: Alternative quantum circuit implementing the dis- 
crete variable DJ Algorithm. Note the absence of the target 
qubit and the use of the operator U / defined in Eq |2.5l 



rithm is usually studied via its corresponding quantum 
circuit. We present a standard circuit version [Tt| in 
Fig.ffl 

The state represented by the lower line in Fig. [1] is re- 
ferred to as the target qubit. In order for easier adaption 
of this circuit to the CV setting, we choose an alternative, 
and equivalent, circuit formulation — one without the 
target state. We take this approach to avoid some of the 
difficulties the target state introduces in TS]. The uni- 
tary operator associated with the oracle function changes 
slightly in this alternative circuit. We discuss these dif- 
ferences before proceeding with analysis of the circuit. 

This simpler algorithm without the target qubit is 
given in Fig. [2j Oracle application is the critical part 
of the algorithm. The oracle construct originally pro- 
posed by DJ is expressed, for x S {0, 1}" and y £ {0, 1}, 
as 



U f :\x)\y)^\x)\y®f(x)). 



(2.3) 



This construction yields a matrix representation for 
the Uf as a permutation matrix, hence always unitary [l[ . 
With respect to the ordered basis 

B = {|0...0)|0),|0...0)|I),...,|I...I)|0),|I..-I)|I)}, 

the unitary matrix U / can be expressed in the following 
insightful form 



Uf 



X f(°-v 



\ 














... x/(' x - 



(2.4) 



with X the 2x2 NOT operator in this case. Here Uf is 
a 2("+!) x 2^" +1 ' matrix, which results from there being 
2" strings (the arguments of /) and an additional target 
qubit. 

The operator Uf can also be expressed in the alterna- 
tive ordered basis 

S' = {|0---0)|-),|0---1)|-),...,|1...0)|+),|1---1)|+)}, 
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We express this as 



Uf 



U f 
1 



with 1 the 2™ x 2™ identity operator. Furthermore the 
operator Uf is expressed as the 2™ x 2™ matrix 



Uf 



/ (_l)/(o-o) o 

(_i)/(o-i) ... 



V 



o 







and thus provides a reduced representation for Uf. It is 
apparent that the operator U / acts on a 2™ x 2™ subspace 
of Uf since 

U f =(Uf®\-)(-\)®(l® n ®\+)(+\). 

We make the assumption that if we have the oracle Uf, we 
we also have the oracle U / . We thus conclude that the 
construction employing both control and target qubits 
is not strictly necessary. That is, one could construct 
this algorithm employing the n-qubit control state only. 
Apparently the choice of representation simply depends 
on the nature of the actual physical implementation. 

We now present a step-by-step analysis of the alterna- 
tive circuit presented in Fig. [2] We shall analyze the CV 
circuit in the same steps for cross reference and compar- 
ison. 



1. State preparation 

We use the hat notation in order to emphasize 
that this analysis is of the algorithm presented in Fig. [21 
which employs n-qubit states and not of that presented in 
Fig. [TJ which employs (n+ l)-qubit states . The n-qubit 
input state of the circuit in Fig. [2] is a string of qubits 
prepared in |0 • • ■ 0). The next step in state preparation 
is to place the state |#o) m to a- n equal superposition of 
all computational basis states 

#®"|*o) >-» =2- n/2 \ x ) 

xe{0,l} n 

for H the single qubit Hadamard operator. 



2. Oracle application 

Given the definition of the reduced operator Uf de- 
fined in Eq. 12. 5\ its effect on the equal superposition of 
basis states expressed in the state is to effectively 
encode the TV-bit string z unitarily into the state 



J7/|*i) ' * |# 2 ) - 2-"/ 2 



(_l)/(o-i) 1 



V(-i)/(i-Dy 



(2.6) 



which is a convenient representation. We shall show that 
this representation naturally extends to the CV setting. 



3. Measurement 

Measurement proceeds by first undoing the superposi- 
tion created during the state preparation step. This is 
achieved through the application of the operator U$ = 
H® n , which modifies the state after oracle application 



*7 3 |* 2 > 
The resultant state is 



a;e{0,l}" 



(2.7) 



(2.8) 



We rewrite Eq. (|2.8p with the operator H® n expressed in 
terms of a recursive definition as follows 



|A \ _o-(n+i)/a (H^-V H^) 

1*3/ - 2 lff®("-l) _Jf<8(n-i; 



Given 



/(_l)/(0-0)N 

(_l)/(o-i) 1 



V(-i)/a-Dy 

(2.9) 



(2.10) 



the combination of Eq. fl2.9p and Eq. (|2.10[) allows us to 
see that all of the rows (and columns) of the operator 
if®" have an equal number of positive and negative ones 
except for the first row, which consists entirely of plus 
ones. It is this feature that permits the constant and 
balanced functions to be distinguished in a single mea- 
surement. 

For the two constant cases, Eq. (|2.9p may be expressed 

as 



A i • 

i -i ••• -i 



*3C) = ±: 



"2" 



v 1 ■■■ l ) 



= ± 



W 



(2.11) 



as only the first row does not result in amplitude cancel- 
lation of the 2" constant amplitude components of the 
state |\l/2) ■ Each of the balanced functions result in the 
amplitudes of the state \^2) having an equal number of 
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positive and negative ones. This feature coupled with ac- 
tion of the operator H ® n results in the first component of 
the state \^>z) having zero amplitude for all the balanced 
functions. We express this result as 



|*3B> =± 



( \ 



w 



(2.12) 



where we use the symbol x to represent that the non- 
zero value(s) will land on the other N — 1 components 
depending on which of the (^ 2 ) balanced functions the 
oracle is set to. It is interesting to note that the number 
of rows in the state \^3b) potentially having a non-zero 
value is N— 1 whereas the number of balanced functions is 
exponential in N. This means that many of the balanced 
states can be expressed as real-valued mixtures of the 
computational basis states with the condition that the 
amplitude of the first component is always zero. 

For the final measurement step, we employ the projec- 
tion operator [l[ defined for m G {0, 1}^ as follows 



M„ 



\m){m\ 



(2.13) 



We are only concerned with the first component as dis- 
cussed above, so for the constant cases we have 

Pr[m=(0---0)] = (* 3C |m (0 ... 0) |# 3 c) = 1, (2.14) 

and for all balanced cases we have 



Pr[m= (0- ••())] = 



':SB 



M, 



(0-0) 



*3b)=0 (2.15) 



as required. 

We have completed the study of the quantum DJ al- 
gorithm in a form that allows us to adapt readily to the 
CV setting. Our strategy will be to construct a CV al- 
gorithm analogous to that shown in Fig. [2] and whose 
operator representation is given by 



1*3 



H® n U f H® n \Q---0). 



(2.16) 



This approach is simpler, and we can worry about 
whether or not an implementation will require target 
states when a particular implementation is considered. 
Before delving into the CV algorithm, we present some 
background CV information. 



E. CV Background 



between two position states leads to the Dirac relation 
(x'\x) = S(x — x'), which must be treated carefully. For n 
the size of Problem Q] the target-less quantum D J algo- 
rithm requires n q ubits, which requires a Hilbert space 
of size V = 2" [llil. The Hilbert space for CV prob- 
lems seems quite generous in this respect as it is infinite- 
dimensional. 

In fact the CV Hilbert space is congruent to the space 
of square-integrable complex functions over the real field 
£ 2 (R) [24]. A function / : [a, b] C, for [a, 6] C K, is in 
C 2 (R) if 



\f(x)\ 2 dx < 00. 



The inner product of two functions /, /' is 



ifV) = / f'*(x)f(x)dx, 



(2.17) 



(2.18) 



with positive definite norm and distance metric defined 
by 



= VW>> d(/,/') = ||/-/'||, 



(2.19) 



respectively. 

Typically in CV quantum information discourse, the 
position states \x) are introduced as a basis set of the 
Hilbert space with each \x) an eigenstate of a position 
operator x, with i£l. Unfortunately the state \x) does 
not exist in the Hilbert space; this problem is evident in 
the standard inner product 



(x\x') =5{x-x'). 



(2.20) 



As 5 is not a proper function, position states are not 
proper states. Fortunately the position states are correct 
as a representation; for example f(x) — (x\f) is the posi- 
tion representation of test function / within the context 
of the rigged Hilbert space. Also Eq. (|2.20[) is meaningful 
in the context of distribution theory. 

A rigged Hilbert space is a pair (Ti, $) such that Ti. is a 
Hilbert space and $ is a vector space that is included by 
a continuous mapping into Ti: <I> C Ti. Elements of $ are 
referred to as 'test functions', and the dual to <& is <&* D 
Ti*, for Ti* dual to Ti and $* comprising generalized 
functions, or 'distributions'. The inner product (/']/) is 
in [0, 1] for any /' G $* and for any f G $ [2ljj . 

Note that the adaptation of the DV DJ algorithm to 
the CV regime needs to be done in the context of a 
computational problem. Here the relevant problem is 
still Problem [TJ and the notion of the oracle remains 
unchanged. Thus, in the CV case, our task is still to 
determine whether the function f z belongs to the set of 
constant functions or to the set of balanced functions. 



The transition from DV to CV quantum information 
requires an extension of Hilbert spaces to rigged Hilbert 
which allows the use of position states 
but restricts dual states to so-called 'test 
An inner product between position states 



spaces 
with x 



21 



G 

functions'. 

and test functions is meaningful, but the inner product 



III. CV REPRESENTATION OF THE DJ 
PROBLEM 

We begin by giving a strategic overview in order to 
convey the key concepts of our approach to developing a 
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FIG. 3: Illustration of the concept for encoding an iV-bit 
string in a region of momentum extending from — P to +P 
using the N = 4, z — 0101 example. Note that each of the 
bits Zj are uniquely represented. 



CV computation model. We follow this with a subsection 
giving some preliminary definitions allowing us to set the 
stage for detailed analysis. We then proceed with a step- 
by-step analysis of our CV DJ algorithm. 
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FIG. 4: An illustrative overview of the four stages of our 
conceptual CV DJ algorithm: (a) The probability distribu- 
tion of the input state wave function positioned at x = xq 
is that of a sine function, (b) The Fourier transform of the 
position state wave function is a "pulse" function in the mo- 
mentum domain, which acts as the encoding "substrate" . (c) 
The iV-bit string z = • ■ • 01 • ■ • 1 modulates this momentum 
"substrate" . (d) The inverse Fourier transform of the encoded 
"square wave" produces a "generalized" sine function whose 
infinite position extent necessitates an optimal measurement 
"window" parameterized by ±<5. 



A. Strategy Overview 

Although we are now working with CV, instead of DV, 
quantum information, the computational problem to be 
solved remains Problem [TJ In other words, we want to 
learn whether the function f z is constant or balanced 
with as few oracle queries as possible. Another way to 
think of this is that we wish to determine the index b G 
{0, 1} such that f z G St,. We now give a conceptual 
overview of our model for CV quantum computation of 
the DJ problem, which we follow later with a rigourous 
treatment. 

In our model of CV quantum computation, we will use 
the continuous position and momentum variables of a 
particle. For R, we use the particle's position wave 

function, 4>(x), to describe where the particle is concen- 
trated and the particle's momentum wave function, </>(p), 
to describe its momentum distribution. The position and 
momentum wave functions are Fourier transform pairs, 
and the relationship between the particle's position and 
its momentum is governed by Heisenberg's uncertainty 
principle. 

There are many position and momentum wave function 
pairs on which we could base our computational model. 
We select our particular pair as follows. First, we wish to 
encode the unknown Af-bit string, z, in the momentum 
domain. We do so because encoding in the momentum 
domain is the continuous analogue of the discrete case, 
where encoding is performed on an equal superposition 
of computational basis states. Second in order to fix one 
of the degrees of freedom of the problem, we want each 
of the bits comprising the string z to be unambiguously 
represented in the momentum space. By unambiguous 



we mean that each of the bits are represented by equal- 
sized, non-overlapping, contiguous regions in the momen- 
tum space. 

Since we want each of the Ambits comprising the string 
to be represented unambiguously, we naturally think of 
each bit as being manifested by a finite-width square 
pulse whose position in momentum space represents the 
bit position in the string z and whose magnitude rep- 
resents the bit value. Continuing along this line of rea- 
soning to the representation of the entire string z, we 
can imagine we have a region of momentum extending 
from — P to +P. All the contiguous momentum pulses 
within this region thus have "width" S p — 2P/N, and for 
j G {0, A — 1}, the j th momentum pulse is centred at 
position — P+ (j + l/2)5 p and takes on value (— l) Zj . We 
illustrate this concept in Fig. [3] for a particular N = 4 
case. 

The picture that thus emerges is that each of the 2 N 
possible strings may be represented by a uniquely shaped 
"square wave" having extent ±P comprising A pulses 
each of width 2P/N and having magnitude ±1. With 
the encoding concept clear, we conceptually illustrate the 
four key stages of the algorithm in Fig.[U We begin with 
a position wave function centered at x = xq and illus- 
trated in Fig. [4] (a). Note that this position wave function 
is a sine function since sinc/pulse functions are Fourier 
transform pairs. In Fig. [¥] (b), we present the momen- 
tum wave function, a pulse function, which acts as the 
"substrate" into which the A-bit strings are encoded. In 
Fig. 0] (c) , the pulse function is encoded with the partic- 
ular A^-bit string z — • • • 01 • • • 1. Finally, the inverse 
Fourier transform of this "square wave" is presented in 
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Fig. [4] (d) . Since the inverse Fourier transforms of finite 
pulses in the momentum domain have infinite extent in 
the position domain, we need to limit the extent of our 
measurement to ±6. 

In summary, we see that our algorithm will need the 
parameters N, P and 5. We note that as N gets large, 
the individual pulse width associated with a single bit 
becomes small appearing to pose a limit on the maximum 
value of N. We will return to this issue once we have 
determined the relationship between P and 8. 

There are many potential models for quantum com- 
putation in a CV setting. We have chosen to study one 
where we unambiguously encode an N-bit string into the 
continuous momentum variable of a particle. Within the 
constraints of this model, we will show that the CV DJ 
is necessarily probabilistic and prove an upper bound on 
the query complexity of the CV DJ problem. 

We speculate that we can't do better than this. For ex- 
ample if the momentum/position pair are described by 
Gaussian/Gaussian functions, as would be the case for 
the physically meaningful states of quantum optics, im- 
perfect encoding of the iV-bit string in the momentum 
domain will result in increased position error. Whether 
or not this will in turn impact the "big Oh" representa- 
tion of the query complexity requires further research as 
does a general proof of a lower bound. The challenge will 
be to show that another strategy can do better than the 
model described herein. 



B. Algorithm Preliminaries 

We now proceed to formalize some of the concepts pre- 
sented in the previous subsection. Here we describe a 
'natural' way of encoding a finite-dimension, iV-bit string 
in a continuous domain. We define the following function, 
along with its Fourier dual, to help us achieve this end. 

For P > 0, the 'top hat' function 

n(p;P,P ) = ( P \n(P,P )} 

1 f l, if P£ [Po-P,Po + P] 
^2P\ 0, if pi [P -P,P +P] 

(3.1) 

will be especially useful in bridging the gap between DV 
and CV quantum information because 

hmn( P ;P,P )=S(p-P ), (3.2) 

so the state |n) is, in some sense, a momentum eigen- 
state \p = Pq) in the limit P — ► 0. The inverse Fourier 
transform of the function e lx ° n (p; P, Pq) is 

4>{x) ={x\4>) 

=7r" 1/2 sinc(P(x- x );P ) 
_ e' p »sin (P(x-xp)) 
yPn(x — Xq) 



where xq defines the position of the sine function. The 
limit oicj)(x) as P goes to oo yields 5(x— xq). The position 
eigenstate \x = xq) is likewise formed in the limit P — > 
oo. 

Now imagine we want to sum a contiguous string of 
"pulses" described by the top hat function (|3.ip with all 
pulses having width S p and the j th pulse having complex 
amplitude ipj. This results in the composite function 

i>(p) = n (p; - p + A< - p + (J + ■ 

3 

This function can also be used as a basis [24| of CV kets 
in Dirac notation as 

/oo 
dp(p\i>)\p), 
-oo 

thus allowing us to encode quantum information in the 
CV domain. Note that ip(p) is the complex amplitude for 
real-valued p. This affords a consistent way of encoding 
a discrete wave function over a continuous domain. 

Before proceeding with a formal analysis of the algo- 
rithm, we give an overview of our proof strategy. The 
oracle is either set to one of two constant strings or to 
one of (jy/ 2 ) balanced strings. A string and its comple- 
ment have indistinguishable probability distributions, so 
there are a total of one constant probability distribution 
plus ^ (jijy 2 ) balanced probability distributions represent- 
ing the possible oracle settings. In order to simplify the 
analysis, we wish to replace this exponential number of 
balanced probability distributions with a single "worst- 
case" balanced probability distribution. Thus we seek a 
particular balanced string (and its complement) whose 
probability distribution is most likely to "fool" us into 
concluding it is a constant string. 

Intuitively, the balanced strings that have the fewest 
number of changes between adjacent bits in the interval 
[— P, P] will be the most "constant like" of the balanced 
strings. There are no balanced strings with zero changes - 
this is the key feature that separates the constant strings 
from the balanced strings. There is however, a single 
pair of balanced strings having only one change. These 
strings exhibit the feature that the first N/2 bits are con- 
stant and the second N/2 bits are the complement of the 
first. We call these strings the anti-symmetric balanced 
(ASB) strings. One of these two strings is illustrated in 
Fig. H](c). Note that all other balanced strings have more 
than one change. 

Our proof strategy begins by making the assumption 
that the ASB case is the "worst case" of all balanced 
cases. We use this assumption to determine the optimum 
value of S, which is the extent of our measurement in 
the position domain and is illustrated conceptually in 
Fig. 2] (d). Given this optimum value of <5, we then prove 
by induction that the worse balanced case is indeed the 
ASB case. 
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FIG. 5: Quantum circuit implementing the CV DJ algorithm 
without the use of the target state. 



1. State preparation 

We have argued previously that we need the Fourier 
transform of the input state to be the top hat function 
defined in Eq. (|3.1|) . We add several conditions that do 
not take away from the generality of the solution. First, 
we want the top hat to have zero phase, which gives Xq = 
and to be centred at Po = 0. Second, we want the pulse 
to have extent ±P. This gives us the simplest form of 
the sine function for the initial state 



(f> (x) = 



sin(Px) 



(3.5) 



C. The CV Quantum DJ Algorithm 

Our strategy is to create a CV analogue of the alterna- 
tive formulation of the discrete D J algorithm presented in 
Fig. O The CV extension of this is presented in Fig. [5] 
Our construction of a CV DJ algorithm employs some 
of Braunstein and Pati's techniques [l8( and avoids the 
pitfalls. In particular, we employ position states as a 
logical representation (states in $*) analogous to the dis- 
crete computational basis states. Encoding is not, how- 
ever, into the position states but rather into test func- 
tions f z € $ with z £ {0,1}^. Furthermore we employ a 
Fourier transform to operate as a CV version of the DV 
Hadamard transform (extending the Hadamard transfor- 
mation to the CV case is not unique [2|. |25|). 

For x the canonical position and p the canonical mo- 
mentum, the Fourier transform maps a function 4>{x) to 
its dual <j){p) according to [26| 



F : <j)(x) i-> <t>(p), 



such that 



and 



Hp) 



4>{x) 



dxe ipx cj>(x), 



dpe- ipx 4>(p). 



Note that we make use of the momentum variable p as 
the Fourier dual of the position variable x. The func- 
tion <f> can be a test function in $ and 4>(x) is the inner 
product of <f> with the position in $*: <j)(x) — (x\cf>). The 
momentum state \p) is the Fourier transform of |x) and 

m = w)- 

With these concepts in order, we now proceed through 
the CV DJ algorithm analogous to the three steps in the 
DV DJ algorithm. The function notation 4>{x) and cf>(p) 
is more convenient here rather than the Dirac notation 
in the previous section. 



We note that the limit of <fio(x) as P — > oo gives a S(x). 
Thus we can think of the quantity P as playing the role 
of the standard deviation in a Gaussian distribution. 

The final step in state preparation is to perform the 
Fourier transform, which yields the top hat function with 
extent ±P 



i 



i, if P e[-P,P] 

2P I 0, ]£p$[-P,P]. 



(3.6) 



This function forms the raw substrate, which will be 
'modulated' by the individual iV-bit strings z. 



2. Oracle application 

We perform encoding by partitioning the real num- 
bers representing momentum into non-overlapping, con- 
tiguous and equal-sized bins. In this digital-to-analogue 
strategy, the width of each p-bin is 2P/N, and 



(3.4) nf } (p) 



W-l-O 



0. otherwise. 



'1-2 



N-i 
N 



The oracle encodes the index z into the function f z as 
follows: 



N-l 



/r ) ( P ) = E(- 1 ) Zl nT(p) 



(3.8) 



where the factor (— l) Zi serves to modulate the phase of 
the top hat function according to the bit value. 

Example 1. Consider the case n — 2; hence N = 2 2 = 
4. As one case, the function corresponding to the four-bit 
string 0011 is 



rf°(p)-nr(p). (3-9) 



The only two four-bit strings yielding constant functions 
would be 0000, for which the function is identically unity 
over the whole domain [—P,P], and 1111, for which the 
function is identically —1 over [—P,P]. Four cases are 
presented in Fig. We refer to the function /ooii(p) o,s 
the "lowest- order" antisymmetric balanced wave as it has 
just one zero crossing in [—P,P]. 




In the limit that N — > oo with P fixed, i— > 
8{p ~ for the midpoint of the i th bin. The limit 
N — > oo thus gives a prescription for approaching a con- 
tinuous variable representation where the z index seems 
to approach a continuum; however this limit yields a 
countable, rather than uncountable, set {z}, and the fi- 
nite domain [— P, P] has important ramifications on the 
nature of the functions corresponding to Fourier trans- 
forms of \li(p). We express the state after encoding as 



(3.10) 



where we observe the "modulating" effect of the encoded 
string f z on the momentum "substrate" <fi(p). 

In the context of the digital-to-analogue strategy, the 
constant functions are analogous to direct current (DC) 
signals and the balanced functions to alternating current 
(AC) signals. The number of zero-crossings corresponds 
to frequency information, and the question of whether 
the output is balanced or constant is essentially a prob- 
lem of querying whether there is a non-zero frequency 
component of the output signal. As noted previously, 
the ASB function has the lowest frequency component. 
We now proceed to analyze the measurement stage. 



3. Measurement 

We have the strings z £ {0, 1}^ encoded into the mo- 
mentum state (|3.10l) . The next step prior to the final 
measurement is to take the inverse Fourier transform of 
this pulse train. For Zj the j th bit of z, this is expressed 



N 



x (e^)^_e*e 



•-)Px 



(3.11) 



FIG. 7: The phasors e* v i x defined in Eq. (pTT2jl for N = 8 
(a) x = 7r/2 phasors range between ±77r/16 in steps of 7r/8, 
and(b) x = 7r/4 phasors range between ±77r/32 in steps of 
it/16. 



The expression given in Eq. (|3.1ip can be simplified to 
yield 



sin(Px/N) 




where we have defined 



n - (2j - 1; 

N 



Px. 



(3.12) 



(3.13) 



We see that the magnitude of an individual generalized 
sine function, (fii^ (x) , is determined by a vector sum 
of N phasors, which is modulated by a particular A^-bit 
string z. 

Note that the phasors, e lip ^ x \ are equiangular divi- 
sions of the angular interval 

[+(N - l)Px/N, —(N - l)Px/N], 

and they exhibit the pairwise complex conjugate property 
<Pj(x) = —(fiN+i-j(x). In Fig. [71 we present the phasors 
for N — 8 with x — ir/2 and x — ir/4 to illustrate these 
features. Note that the phasors are added constructively 
or de-constructively depending on the phase of the an- 
gles, which results from the term (— l) Zj . This effect 
defines the magnitude of the resulting sine function. 

We note that the only functions with (f>z (0) ^ 
are the two constant sine functions. This is clear given 
that ^2^ = i{—l) Zj = for the two constant cases, and 
Ylf=ii~^) Zi = since for all balanced cases, the latter 
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sum always resolves to N/2— N/2 = 0. This feature of the 
set of (j^ 2 ) + 2 sine functions represented by <pi^ (x) de- 
fined by Eq. (13.12j) implies the strategy for measurement 
that will distinguish between the constant and balanced 
cases. 

In order to refine this strategy, we focus on two cases. 
The first of these cases is for the two constant functions 
for which Eq. (|3.12p gives the probability distribution 



0.35 r 



0.35 r 



V c {x) = \<^\x)\ 2 = 



sin 2 (Pa;) 
Pttx 2 '' 



(3.14) 



where we have 







C € 



The second case deals with the two balanced functions 
having the lowest 'frequency' content, which occurs when 
the first N/2 bits and the last N/2 bits have opposite 
values. 

We think of these two balanced functions as having 
the lowest frequency content since Eq. (|3.8[) has a single 
zero crossing in the interval [— P, P] for these two bal- 
anced strings only. All other balanced strings have more 
than one zero crossing and thus higher frequency con- 
tent. For this pair of balanced functions, which we call 
the antisymmetric balanced (ASB) functions, Eq. (|3.12|) 
gives the probability distribution 



V^{x) = \4>lll{x)\ 

where we have 



W/_\i2 _ (cos(Pa;) - l) 2 



Pttx 2 



ASB e < - • - 01 — 1, o_^_o 

N/2 N/2 N/2 N/2 



(3.15) 



(3.16) 



Note that of the {jy 2 ) balanced functions, there are many 
that are also antisymmetric about the midpoint. How- 
ever, we reserve the term ASB for these two lowest-order 
antisymmetric balanced functions. 

We will use these cases to bound the success proba- 
bility of distinguishing between the constant and all bal- 
anced cases. We first illustrate the concept of frequency 
in the following example. 

Example 2. Again consider the case P — 1, n = 2; 
hence N = 2 2 = 4. As one case, the function corre- 
sponding to the four-bit string 0011 is 



u oon 



,(4) = Bfo(E/4) 

—i (coscc 



e 4 + e 4 



(3.17) 



This function corresponds to the N = 4 ASB function. 
The probability distributions for the four distinct N 
1 cases are presented in Fig. [31 We clearly see that of 
the three balanced cases, the N — 4 ASB function has 
probability peaks closest to Xq = 0. 



FIG. 8: The probability distributions |<?i>i 4) (x)| 2 for (a) z = 
0000, (b) z = 0011, (c) z = 0101, and (d) z = 0110. We 
clearly see that of the three balanced cases (b) through (d), 
the ASB function (b) has probability peaks closest to xq — 0. 



Our measurement strategy is to measure the proba- 
bility distribution in a small band around the position 
xq = parameterized by ±5. The CV analog of the 
projection operator given in Eq. (|2.13p is defined as 



where 



E. 



D b a (x) 



D b a (x)\x)(x\dx, 



1, if a < x < b 
0, otherwise. 



(3.18) 



(3.19) 



Due to the symmetry of the sine functions about xq , we 
set a — —S and b = +5. We now need to determine the 
optimal value of S that will maximize our ability to distin- 
guish between the constant and balanced cases. We will 
determine the optimum value of S by first assuming that 
the probability distribution Vasb (x) given by Eq. (|3.15|) 
dominates all other balanced probability distributions in 
the region [—5,5]. After using this assumption to deter- 
mine a value for the optimal delta, we will state and prove 
a theorem justifying our assumption. As an illustration 
that our assumption is true for the N = A case, we plot 
the four distinct cases in Fig. [HI 

The ability to effectively distinguish between two ran- 
dom events is proportional to the separation of the indi- 
vidual probabilities of occurrence. Thus we need to select 
5 such that we get as much separation between the con- 
stant distribution and the ASB distribution as possible. 
Given this concept we can think that when we make a 
measurement we are distinguishing between two events, 
the probabilities for which we define as follows 



Pr ConBt (5)=Pr 

and 
Pr ASB (5) =Pr 



4 N) 



dN) 



Vo{x) 



Vasb{x) 



= E s _ s {V {x)) 



(3.20) 



E s _ s (V ASB (x)). 

(3.21) 
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FIG. 9: For P = 1, the optimal value of 5 = f . This 
graph shows that only the Constant and the Antisymmet- 
ric Balanced Functions significantly contribute to probability 
between ±5. 



We can determine the optimum value of S by maximiz- 
ing the expression |Pr Const ((5) — Pr ASB (<5)|. It suffices to 
find the value of S for which ^|Pr Conat (5) — Pr ASB (<5)| = 0, 
which may be expressed as 



FIG. 10: Definition of the phasor angles for the N = 4 base 
base. Note that the effect of varying x over [— tt/2, tt/2] simply 
focuses or expands the double angles 2ipi and 2i^2 proportion- 
ally. 



^|Prco n3t (<5)-Pr AS B(<5)| 



and 



sin 2 (Pa;) (cos(Px) - if 



g \ Pttx 2 Pttx 2 
sin 2 (PS) (cos(PS) - l) 2 



d.r 



PnS 2 



PnS 2 



0. 



(3.22) 



This occurs where cos(P<5) = cos(PS) 2 for 5^0, which 
gives a global maximum at 8 = It is interesting to 
think of this result as an uncertainty relationship 



PS = 



(3.23) 



We shall return to this concept in our discussion in the 
conclusion. We have determined the optimum value for 
S based on our assumption that for — S < x < S the 
balanced probability distribution V AS13 (x) dominates all 
other balanced probability distributions. We now pro- 
ceed to prove this assumption. 

In order to proceed with the proof, we define a set 
<f> of m pairwise conjugate angles with 2m = N. Note 
that N is not restricted to being equal to 2 n for the pur- 
pose of this proof. Also for the purpose of this proof, 
we set P = 1 and incorporate x into the definition of 

<Pj = ( N ~ { n~ 1) ) x for _7r / 2 - x - ^l 2 - Wc lct $ = 
{ip!,(p 2 , ■ ■ .,<p m > < Pm+l) ■ ■ ■>ftm} with j = 1 , . . . , in and 
note the pairwise conjugate property ifj — — </?2m+i-j- 
Now consider S = Sj=o 9{j) e%Vi where g : [2m] i— > ±1 
subject to the balanced condition J2j = 0, then 

Theorem 1. Max |S| occurs under the specific balanced 
conditions 



g(J) = 



1 if ' 1 < j < m 
— 1 if m + 1 < j < 2m, 



g{j) 



— 1 if 1 < j < m 
1 if m + 1 < j < 2m, 



(3.25) 



which we refer to as the asymmetric balanced functions 
(ASB). 

Proof. Proof is done by induction on m. We begin with 
the base case m— 1, N = 2. This case is trivial since the 
only balanced cases are the two ASB cases represented 
by the strings {01, 10}. We proceed with the base case 
for m = 2, N = 4. This case is a little more involved. 
We begin by labelling the angles and phasors as shown 
in Fig. [TO] 

There are (2) = 6 balanced cases, and we have to 
consider the strings {0011,0101,0110,1100,1010,1001}. 
Since the latter three are complements of the first three, 
we have to consider only three vector sums. 

With reference to Fig. [TOl we have S{ on} = e lipi + 
e lV2 — e~ tv>1 — e~ lV2 . We simplify and express the resultant 
along with the three other cases as 

1. Si = S { {oou},{iioo}} = ±2«(sin(</Ji) + sin(y> 2 )) 

2. S 2 = S {{0101}j{1010}} = ±2(cos(<y3i) - cos(v? 2 )) 

3. S 3 = S {{ oiio},{iooi}} = ±2«(sin((/3i) - sin(y> 2 ))- 
Clearly |Si| > |S 3 |. We use the trigonometric identities, 



|Si| = 2sm ( I cos I 

« I .) • [ >1 + ( PA . (<Pl - r'2 

I S 2 1 = 2 sin I I sin 



(3.26) 



(3.24) to establish the relationship between |Si| and |S 2 |. We 



note that max ( ipi 2 y2 



jj) x and mir 
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x for all N and the specified range of x. Since cos sc > 
sinx for 0<a;<l/2, we conclude that 



^2 \ ^ . ( <fl 

> sin 



<P2 



and thus | S x | > | S 2 1 for < x < n/2. This proves that 
the theorem is true for the m = 2, N = 4 base case. We 
are now ready to prove the inductive step. 

We consider two cases. Case (i) assumes every pair is 
balanced. By this we mean that g(j) — —g(2m + 1 — j). 
By inspection, this gives the same result as | S x | and |S 3 | 
for the in = 2, N = 4 case. Case(ii) assumes that Case(i) 
is not true and is proved by induction. Since Case(i) 
is not true, there must exist two non-balanced pairs for 
which g(j) — g(2m + l—j) = +1 and g(k) = g(2m + 1 — 
k) = — 1. As an illustration in the m — 4, N = 8 case, 
the balanced string {01000111} has this property. The 
inductive step is 



$>C7)« 



3 = 1 



< |S ({/, 2m + 1 - I, k, 2m + 1 - k}) \ 
+ \S(w)\, (3.27) 



where S(w) is maximized for the m = 2, N = 4 base 
case. Only when \S ({I, 2m + 1 - I, k, 2m + 1 - k}) | it- 
self is maximized is equality achieved and the total sum 
maximized. This occurs for the ASB strings. □ 

We have established that we can bound the probabili- 
ties of determining whether an unknown function is bal- 
anced or constant in a single query in a CV setting. In 
the next section will determine an upper bound for the 
query complexity of a CV algorithm in terms of success 
probability in terms of the number of queries. 



IV. BOUNDING THE QUERY COMPLEXITY 
OF THE CONTINUOUS VARIABLE DJ 
ALGORITHM 

Before we bound the query complexity, we make some 
important observations regarding the comparison be- 
tween the discrete DJ algorithm and the CV DJ algo- 
rithm. 

First, we note that probability distributions V c {x) and 
"Pases (x) defined by Eqs. (I3.14|) and ()3.15|) respectively, 
are in Ti.2, the Hilbert space of £ 2 (R) functions over the 
interval [—00,00]. This implies that since we are measur- 
ing over a finite interval, the CV DJ algorithm is neces- 
sarily probabilistic. Furthermore, we noted that P and 8 
are related by the uncertainty relation given in Eq. (|3.23D . 
This leads to the conclusion that even in the limit of the 
improper delta function 8 (x — xq), the CV DJ algorithm 
remains probabilistic. This conclusion is contrary to that 
made in [3j . 



Second, we compare the operator descriptions of the 
DV DJ and the CV DJ, which we express as 



|*3) = H® n U f H® n \0 •■■0) 



(4.1) 



The first equation represents the quantum DJ algorithm 
operator expression given in Eq. (|2.16|) . The second equa- 
tion is the analogous CV operator expression determined 
by concatenating the steps of the previous section. There 
is a high degree of similarity between these two expres- 
sions, but there are mathematical subtleties. 

1. The CV position state 4>o(x) is not a perfect analog 
to the computational basis state |0 ■ • • 0) except in 
the limit. However, this limit creates a state that 
is not in the RHS we argued is necessary for con- 
sistency (2lj . 

2. The continuous Fourier transform is not equal to 
the CV extension of the Hadamard operator in 
a CV parameterized system with a finite Hilbert 
space. It is however, a convenient extension when 
the Hilbert space is infinite. 

3. Finally, the diagonal operator U* given by Eq. (|2.5p 
has each entry taking on the value ±1 dependent 
on the value of f z . The CV analogue to this oper- 
ator is the function fz' (p) , where each of the N 
partitions of the real interval [— P, P] similarly take 
on the value ±1. 

We now determine numerical values of the probabilities 
determined in Eqs. (|3.20p and (|3.2ip . We can readily 
calculate the probability of detecting if the function is 
constant 



PlV 



5 sm 2 {Px) 



dx 



_ cos(2<SP) + 2SPSi{28P) - 1 

where the sine integral is given by 

Si(z) = / At. 

.In t 



(4.2) 



(4.3) 



Note this probability depends only on the product PS. If 
the function is the lowest-order antisymmetric balanced 
(ASB) 

_ f s (cosjPx) 1)2 
FrASB _ J_ s Pirx* dX 

_ -8 sin 4 (^) + ASPSi(SP) - 2 ( 5PSi(2<SP) 



SPtt 



(4.4) 



For PS = 7r/2, the numerical values of these two proba- 
bilities are 



Pr r 



2(?rSi (tt) - 2) 



0.77, 



(4.5) 
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and 

P rASB ^ 47rSi ^- 2 27rSi ^- 4 .0.16. (4.6) 

Given this probabilistic nature of the CV DJ algorithm, 
we need to develop a strategy to bound the error prob- 
ability. We will employ the technique sometimes called 
probability amplification (27l . |28| . 

Our strategy will be to make m repetitions of the CV 
DJ algorithm where we assume that the oracle is set to 
the same function for each of the repetitions. Each rep- 
etition ends with a measurement. From this sequence 
of measurements we want to determine whether the un- 
known function is balanced or constant with high prob- 
ability. 

Theorem 2. An error of 0(e~ m ) can be achieved by 
making O(m) repetitions of the CV DJ algorithm. 

Proof. We will adopt the convention that when we make a 
query to the CV D J algorithm we either detect something 
(algorithm returns a 1), or we do not (algorithm returns 
a 0). We can thus treat multiple queries as a sequence 
of Bernoulli trials . We assume that we have set our 
measurement limits to the optimal ±<5. The two events 
we are trying to uncover are the constant cases where, 
for ease of calculation we set the probability of detecting 
something is Pr c > 3/4, and the balanced cases where the 
probability detecting something is Pr B < 1/4. Note that 
we have set the probabilities to these rational numbers 
for illustrative purposes and to simplify the calculation. 
We can make this arbitrary setting, and we will get the 
same result as long as the probabilities are bound from 
1/2 by a constant. 

If each measurement is based on an independent prepa- 
ration of the state <po(x), then each of the queries are 
independent. After a series of m queries, we can use 
the Chernoff bounds of the binomial distribution to am- 
plify the success probability [28|, Hi| . The simplest (but 
somewhat weak) Chernoff bound on the lower tail is given 
by [H as 

Pr[X < (1 - e)fj,)\ < e~^, (4.7) 
and on the upper tail as 

Pr[X> (1 + e)//)] (4.8) 

where /1 is the expected mean of the resulting binomial 
distributions after m queries, and e is the relative dis- 
tance from the respective means. 

First, we bound the lower tail corresponding to the 
distribution of the constant case for which we have [i = 
mV c - Here we set e = i, which expresses the probability 
for the value being less than half way between the two 
means as Pr[X < (m/2)] < e~24. Similarly, we bound 
the upper tail for the balanced case for which we have /i = 
m V B . Here we set e = 1 , which expresses the probability 
for the value being greater than half way between the two 



means as Pr[X > (m/2)] < e~T5. Clearly the success is 
worse for the lower tail allowing us to bound the success 
probability of the CV D J algorithm after m queries as 

Pr[Success] > 1 - e~^. (4.9) 

This gives an error probability that is O (e~ m ) as re- 
quired. □ 

We note that this is of the same order as the exponen- 
tially good success probability we have for the classical 
probabilistic approach given by Eq. (|2.2[) . Also note that 
this query complexity is independent of the value of N . 
We have made no attempt to obtain a tighter bound pre- 
ferring to show only that the success probability of the 
CV D J algorithm is of the same order of that of the clas- 
sical probabilistic approach to solving the DJ problem. 

V. CONCLUSIONS 

In this paper we have presented a rigorous framework 
for the analysis of the DJ oracle identification problem in 
a CV setting. The rigged Hilbert space (RHS) affords a 
consistent transition from the traditional discrete Hilbert 
space to the CV setting. Our framework allows us to 
define a consistent way of encoding AT-bit strings into 
functions over the real numbers. 

We have used this framework, and the selection of the 
sinc/pulse Fourier transform pair, to prove that a CV 
implementation of the DJ algorithm cannot provide the 
exponential speed-up of its discrete quantum counter- 
part. Additionally, we have presented a bounded-error, 
upper bound on the query complexity of the DJ prob- 
lem within the constraints of our model. The lack of 
speed-up results from an uncertainty principle between 
the ability to encode perfectly in a continuous represen- 
tation and the subsequent inability to measure perfectly 
in the Fourier-dual representation. This uncertainty re- 
lationship is manifest in Eq. (|3.23[) . which relates P, the 
encoding extent, to <5, the measurement extent. A natu- 
ral extension of this work would be prove a lower bound 
perhaps exploring the techniques along the lines of (30j 
from the perspective of different Fourier transform pairs. 

This uncertainty relationship appears to be a natural 
feature of the CV setting, but it could also be used to 
advantage. There is likely to be oracle function symme- 
tries that are particularly suited to different CV settings. 
For example in Sec. IIIIl we showed that balanced func- 
tions with a higher number of zero crossings create sine 
functions with frequency components further away from 
xq. It appears that an oracle identification problem de- 
signed to separate balanced functions according to fre- 
quency separation could be implemented in a CV setting 
and possibly provide advantage over classical or discrete 
quantum settings. 

Furthermore, it would be interesting to classify the bal- 
anced functions from the perspective of different coherent 
states [3l| in CV parameterized settings of both finite and 
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infinite dimensions. The former would naturally involve 
the study of the coherent spin systems [HJ • Furthermore 
the use of squeezed spin states should be studied [33j ] - 
Infinite dimension systems would naturally involve the 
study of implementations involving the coherent states 
of quantum optics [H, [34[ ■ 

Additionally, we have set up this framework in a man- 
ner that should allow any oracle identification problem to 
be analyzed in a similar manner in the CV setting. An 
implementation of a discrete quantum oracle, for exam- 
ple [iH, [3(| , requires a unitary operator representing the 
oracle. Provided we can create a diagonal representation 
of this oracle along the lines of Uf given in Eq. (|2.16p , our 
framework will naturally extend to it. Of course we need 
to be able to create an implementation of these oracles 
and that remains an important open question. 

Other avenues of the extension of this framework in- 
clude CV implementations of other hidden subgroup 
problems. The solution of Simon's problem [371 ] in this 
setting would be an obvious starting point as would 
the exploration of a CV implementation of Shor's algo- 
rithm 38]. Additionally, the CV framework could be 



extended to include analysis of noisy oracles along the 
lines of [11 El- 

In closing we note that the transition from a discrete 
quantum information setting to a CV setting has many 
subtleties. In particular the improper delta functions 
must not be used. Limiting behaviour can be explored, 
but only if the limits are taken from the perspective of 
functions defined in the rigged Hilbert space. 



Acknowledgements 

We appreciate financial support from the Alberta In- 
genuity Fund (AIF), Alberta's Informatics Circle of Re- 
search Excellence (iCORE), Canada's Natural Sciences 
and Engineering Research Council (NSERC), the Cana- 
dian Network Centres of Excellence for Mathematics 
of Information Technology and Complex Systems (MI- 
TACS), and General Dynamics Canada. PH is a Scholar 
and BCS is an Associate of the Canadian Institute for 
Advanced Research (CIFAR). 



[1] M. A. Nielsen and I. L. Chuang, Quantum Computation 
and Quantum Information (Cambridge University Press, 
Cambridge, UK, 2000). 

[2] D. Gottesman, A. Kitaev, and J. Preskill, "Encoding a 
qubit in an oscillator ", Phys. Rev. A 64, 012310 (2001). 

[3] S. L. Braunstein and A. K. Pati, eds. Quantum Infor- 
mation with Continuous Variables, (Kluwer, Dordrecht, 
2003); arXiv: |quant-ph/0207108| 

[4] A. Furusawa, J. L. S0rensen, S. L. Braunstein, C. A. 
Fuchs, H. J. Kimble, and E. S. Polzik, "Unconditional 
quantum teleportation" , Science 23, Vol. 282. no. 5389, 
pp. 706-709, (1998) 

[5] F. Grosshans and P. Grangier, "Continuous variable 
quantum cryptography using coherent states", Phys. 
Rev. Lett. 88, 057902-1-057902-4 (2002). 

[6] J. Appel, E. Figueroa, D. Korystov, M. Lobino, and A. I. 
Lvovsky, "Quantum memory for squeezed light", Phys. 
Rev. Lett. 100, 093602 (2008). 

[7] D. Akamatsu, Y. Yokoi, M. Arikawa, S. Nagatsuka, T. 
Tanimura, A. Furusawa, and M. Kozuma, "Ultraslow 
propagation of squeezed vacuum pulses with electromag- 
netically induced transparency", Phys. Rev. Lett. 99, 
153602 (2007). 

[8] S. L. Braunstein, "Error correction for continuous quan- 
tum variables", Phys. Rev. Lett. 80, 4084 (1998). 
[9] J. Eisert and M. B. Plenio, "Distilling Gaussian states 
with Gaussian operations is impossible", Phys. Rev. Lett. 
89, 137903 (2002). 

[10] S. D. Bartlett and B. C. Sanders, "Efficient classical sim- 
ulation of optical quantum information circuits", Phys. 
Rev. Lett. 89, 207903 (2002). 

[11] S. D. Bartlett, B. C. Sanders, S. L. Braunstein, and 
K. Nemoto, "Efficient classical simulation of continu- 
ous variable quantum information processes", Phys. Rev. 
Lett. 88, 097904 (2002). 



[12] L. Blum, M. Shub, and S. Smale, "On a theory of 
computation and complexity over the real numbers: 
NP-completeness, recursive functions and universal ma- 
chines", Bull. Am. Math. Soc. 21, 1 (1989). 

[13] "Continuous turing machines", Manuscript, not dated. 

(Available at 

http:// www.mathpages.com/home/kmathl35.htm . Ref- 
erence downloaded 03 Dec. 2008.) 

[14] D. Deutsch and R. Jozsa, "Rapid solution of problems 
by quantum computation", Proc. Royal Soc. (Lond.) A 
439, 553 (1992). 

[15] D. Deutsch, "Quantum theory, the Church- Turing princi- 
ple and the universal quantum computer", Proc. of the 
Royal Society (London) A 400, 97 (1985). 

[16] R. Cleve, A. Ekert, C. Macchiavello, and M. Mosca, 
"Quantum algorithms revisited", Proc. Royal Soc. 
(Lond.) A 454, 339-354 (1998). 

[17] D. Bacon, "Lecture 4: Quantum algorithms from sum- 
mer school in Siena, Italy. August 30-September 2, 2005" , 
Manuscript, 2005 .(Available at 
http://www.cs.washington.edu/homes 
/dabacon/teaching/siena. Reference downloaded 20 Oct. 
2008.) 

[18] S. L. Braunstein and A. K. Pati, "Deutsch- Jozsa algo- 
rithm for continuous variables" , in || . 

[19] L. K. Grover, "A fast quantum mechanical algorithm for 
database search", Proc. 28th Ann. ACM Symp. on The- 
ory of Computing (STOC '96), pp. 212-219 (1996). 

[20] A. Ambainis, K. Iwama, A. Kawachi, R. Raymond, and 
S. Yamashita, "Robust quantum algorithms for oracle 
identification", arXiv: quant-ph/0411204 

[21] R. de la Madrid, "The role of the rigged Hilbert space in 
quantum mechanics", Eur. J. Phys. 26, 287 (2005). 

[22] R. Blume-Kohout, C. M. Caves, and I. H. Deutsch, 
"Climbing mount scalable: Physical resource require- 



15 



ments for a scalable quantum computer", Found. Phys. 
32, 1641 (2002). 

[23] C. M. Caves, "Physical resources, entanglement, and 
the power of quantum computation", Lecture given in 
SQuInT Summer Retreat University of Southern Cali- 
fornia, 2005 July 7. Reference downloaded 03 Dec. 2008. 

[24] C. I. Tan, "Notes on Hilbert Space", Manuscript, not 

dated. (Available at 

http:/ /jcbmac. chem.brown.edu/baird/QuantumPDF 
Reference downloaded 20 Oct. 2008.) 

[25] S. D. Bartlett, H. de Guise, and B. C. Sanders, "Quan- 
tum encodings in spin systems and harmonic oscillators" , 
Phys. Rev. A 65, 052316 (2002). 

[26] R. N. Bracewell, The Fourier Transform and Its Appli- 
cations, 2nd Ed. (McGraw-Hill, New York, 1986). 

[27] M. Adcock "The classical and quantum complexity of 
the Goldreich-Levin problem with applications to bit 
commitment" , Master of Science Thesis, Department of 
Computer Science, University of Calgary (2004). 

[28] J. F. Cann y, "Chernoff bounds", Manuscript, 2001 
(Available at http:/ /www.cs.berkeley.edu~jfc/csl74/lecs 
/leclO/leclO.pdf. Reference downloaded 20 Oct. 2008.) 

[29] T. H. Cormen, C. E. Leiserson, and R. L. Rivest, Intro- 
duction to Algorithms, MIT Press, Cambridge, MA 1990 

[30] P. H0yer and R. Spalek , "Lower Bounds on quantum 
query complexity", Bull. Euro. Assoc. for Theoretical 
Computer Science Vol. 87, 2005. 

[31] A. Perelomov, Generalized Coherent States and Their 



Applications (Springer- Verlag, New York, 1986). 
[32] F. T. Arrechi, E. Courtens, R. Gilmore, and H. Thomas, 

"Atomic coherent states in quantum optics" , Phys. Rev. 

A 6, 2211 (1972). 
[33] M. Kitagawa and M. Ueda, "Squeezed spin states", Phys. 

Rev. A 47, 5138 (1993). 
[34] U. Leonhardt, Measuring the Quantum State of Light 

(Cambridge University Press, Cambridge UK, 1997). 
[35] E. Bernstein and U. V. Vazirani, "Quantum complexity 

theory", SIAM J. on Comp. 26, No. 5, pp. 1411-1473 

(1997). 

[36] R. Cleve, W. van Dam, M. Nielsen, and A. Tapp, "Quan- 
tum entanglement and the communication complexity of 
the inner product function" , Lecture Notes in Computer 
Science 1509 (Springer- Verlag), pp. 61-74 (1999). 

[37] D. R. Simon, "On the power of quantum computing", 
SIAM J. on Comp. 26, No. 5, pp. 1474-1483 (1997). 

[38] P. W. Shor, "Polynomial-time algorithms for prime fac- 
torization and discrete logarithms on a quantum com- 
puter", SIAM J. on Comp. 26, No. 5, pp. 1484-1509 
(1997). 

[39] M. Adcock and R. Cleve, "A quantum Goldreich-Levin 
theorem with cryptographic applications" Proc. 19 th In- 
ternational Symp. Theor. Aspects Comp. Sci. (STACS 
2002), H. Alt and A. Ferreira, eds., Lecture Notes in 
Computer Science 2285 (Springer- Verlag), pp. 323-334 
(2002). 



